Meissner-London currents in superconductors with rectangular cross section 
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Exact analytic solutions are presented for the magnetic moment and screening currents in the 
Meissner state of superconductor strips with rectangular cross section in a perpendicular magnetic 
field and/or with transport current. The extension to finite London penetration is achieved by an 
elegant numerical method which works also for disks. The surface current in the specimen corners 
diverges as Z -1//3 where 1 is the distance from the corner. This enhancement reduces the barrier for 
vortex penetration and should increase the nonlinear Meissner effect in d-wave superconductors. 
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The main feature of superconductors is that they expel 
weak magnetic fields H from their interior. This Meissner 
effect was described quantitatively by the London broth- 
ers, who showed that H penetrates exponentially to the 
London depth A Q|, and by Ginzburg and Landau and 
by Pippard, who introduced the superconducting coher- 
ence length £. In extreme type-II superconductors with 
A ^S> £, the correction caused by finite £ to the Meissner- 
London state usually may be disregarded, but now the 
penetration of magnetic flux in form of Abrikosov vor- 
tices has to be considered. Vortex penetration is gov- 
erned by surface barriers which increase with decreasing 
£. Both the microscopic Bean-Livingstom barrier [^]|| 
and the macroscopic geometric barrier | § depend on 
the surface screening currents flowing in the Meissner- 
London state. These currents crucially depend on the 
specimen shape and on the orientation of the applied 
magnetic field H. The screening current is particularly 
large near sharp edges, where it causes a reduction of the 
field of first penetration of vortices and possibly increases 
the nonlinear Meissner effect in d-wave superconductors 
H||. Knowledge of the Meissner-London state of thin 
platelets is also required for a correct evaluation of cer- 
tain precision measurements of the London penetration 
depth A §. 

In spite of its fundamental nature, exact solutions of 
London theory exist only for the trivial (and less impor- 
tant for experiments) longitudinal geometry of infinite 
slabs and cylinders in parallel H, and for the sphere and 
infinite cylinder in arbitrarily oriented H jjj. Even for 
the ideal Meissner state, i.e. the limit A — > 0, the only 
nontrivial solution we know of is the ellipsoid |To 11 and 
strips with elliptic |T|] and oval [jl3| cross sections. In 
general, when A is much smaller than the smallest exten- 
sion of the superconductor, the surface screening current 
is J = \H\\\ where Hn is the component of H(r) at and 
parallel to the surface, and its density is approximately 



j = (J/X)exp(-S/X), 



(1) 



where 6 is the distance from the surface. 

In the present paper we first give an exact solution for 
the surface current and magnetic moment in the ideal 
Meissner state (A = 0) of an infinitely long strip with 
rectangular cross section in a perpendicular magnetic 
field H or with transport current /. We then show how 
the Meissner-London state of this strip (and of circu- 
lar disks) can be computed for arbitrary A. In all these 
cases the solutions for H > 0, / = and H = 0, I > 
may be superimposed linearly to give the general (less 
symmetric) result for simultaneously applied H and /. 
The rectangular cross section of the strip fills the area 
—a < x < a, — b < y < b; the currents flow along the 
strip \\z, and H is applied along y, see inset in Fig. 1. 
The result for arbitrary orientation of H is obtained by 
linear superposition of the solutions with H\\y and H\\x. 

The surface currents of an ideally diamagnetic rectan- 
gular strip in a perpendicular field H\\y or with transport 
current I\\z can be calculated by conformal mapping. We 
give here the main results starting with the case H > 0, 
1 = 0; details will be published elsewhere. Defining the 
function 
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[/ is the sum of two incomplete elliptic integrals, 
f(s,m) = E(6,k) - (1 - k 2 )F{e,k), s = sin0, m = fc 2 , 
0<s<l,0<m<l] we may write the screening cur- 
rents and the magnetic moment in parametric form, with 
curve parameters s and m. First we find m(b/a) from 



b = /(l,m) 
a /(1,1-m) 

and then the currents as functions of x(s) and y(s): 
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The (negative) magnetic moment of the strip per unit 
length along z is 
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(6) 



with m — m(b/a) from Eq. (3). These exact results are 
shown in Figs. 1 and 2. 

Some useful approximations and limiting cases are as 
follows. For all aspect ratios < b/a < oo one has 
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with error |e| < 0.18%. The exact limits are m — > Ah/ita 
(b <C a) and m — > 1 — Aa/itb (b » a). The magnetic 
moment (6) for all ratios b/a is well fitted by 
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with p = 10+| ln(6/a) + 0.288| 1968 and deviation < 
e < 0.004, see Fig. 2. Formula (8) is more accurate than 
the fit given in Ref. 0. The exact limits are (Fig. 2): 
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6< a, 
b > a. 



(9) 
(10) 



The surface currents diverge at the corners symmetri- 
cally as Z -1 / 3 where I is the distance from the corner, e.g. 
I = a — x and I = b — y. Near the corners one has 



J r< 
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(11) 



At the equator J(a, 0) = H / y/m holds and near the poles 
J(x, b)/H « (x/o)/(l, 1— m)/(l— m) which equals nx/Aa 
at b ^> a and i/a at i « n. For long slabs one has 
J(a,y) « H except near the corners [see Eq. (11)]. For 
thin strips (b <C a) Eq. (4) yields 



(12) 



J(x,b) = H x (a 2 - x 2 )~ 1/2 



in agreement with the known sheet current 2 J{x) |15|,|16| 
On the edge of thin strips one has the fit 
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(13) 



with relative error < 1% for y/b < 0.92 or < 2% for 
y/b < 0.97. More limiting expressions and fits are easily 
derived from the exact Eqs. (2-6). 

For the strip with current / > and H = we find 
the Meissner surface currents 
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These expressions are invariant when interchanging a,b 
and x, y, which replaces m by 1 — m. In the corners the 
current again diverges as Z -1 / 3 , cf. Eq. (11): 
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Accounting for a small but finite London depth A <C 
a, 6, the screening currents penetrate exponentially from 
the surface, e.g., j{x, y) = J{x, 6)A _1 exp[(y — b)/X] near 
y = b. The magnetic moment M = Jdx Jdy x j{x,y) is 
thus slightly reduced to approximately 



M(a,b, A) « M(a — A, 6 — A, 0) 



(17) 



with M(a, b, 0) from Eqs. (6, 8-10). Similar formulas are 
valid for specimens of any shape if the radius of cur- 
vature of the surface is much larger than A. For our 
strip, however, the sharp rectangular corners give an ad- 
ditional contribution <5M corner to M(A) — M(0). When 
[A/min(a, 6)] 1/3 < 1, one has from Eq. (11) 



<5M cornor oc {\ 2 /ab) l ' 3 M. 



(18) 



This nonanalytical term dominates in dM/dX and may 
explain some experimental findings in Ref. || . In the op- 
posite case, A>o,ii, the vector potential of the induced 
current density is negligible, thus j ~ —Hx/X 2 and 



M(a, b, A) ps -(4a 3 6/3A 2 )iJ . 



(19) 



Next we show how the current density j(x, y) and mag- 
netic moment M of thick strips (and of disks) for finite 
London depth A can be obtained in an elegant way, avoid- 
ing the calculation and cutoff |l^] of the magnetic field 
around the strip. The static London equation reads 



- AVoj = A = Aj + A q 



(20) 



where Aj is the vector potential of the supercurrent den- 
sity j and A a is the vector potential of the applied field, 
e.g. A a = —ixfioH for strips (z = unit vector along z). 
Inverting /ioj = ~ V 2 A.j we get for thick strips 



^ (rl = _ /(|1 / d V^-^i(r') 



(21) 



with r = (x, y), A — zA, j = zj. The integration is over 
the strip cross section. From Eqs. (20,21) we have 



A a (r)= Mo /dV 
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Solving for j and using A a = —xi-J.qH wc obtain 
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j(r) = H J d 2 r'K(r, r') x , (23) 
^l+A^r-O]" 1 . (24) 



The A dependent integral kernel if(r, r') may be com- 
puted by choosing r and r' on a grid and inverting the 
resulting matrix similar as shown in Refs. |]l4] , p^| . 

From Eq. (23) the current density induced by H in the 
Meissner-London state is obtained by a simple integra- 
tion over x and y. Due to the symmetry of j{x,y) = 
j(x,-y) = -j(-x,y) = -j(-x,-y) it suffices to inte- 
grate over a quarter of the strip cross section, < x < a, 
< y < b, if the kernel is made symmetric. 

Similar equations follow for the strip with current /, 



j(r) = 1 J d 2 r'K{v,v') j J d 2 r J d 2 r'K(r,r') 



(25) 



with the same kernel, Eq. (24), which however has a dif- 
ferent symmetric form since now j(x,y) — j(x,—y) — 
j{—x,y) = j{—x,—y) holds. For a strip with both ap- 
plied field H and current / the two solutions, Eqs. (23,25) 
may be superimposed. The same method works for thick 
disks in axial field if the appropriate kernel K is used, 
similarly as shown in Ref. (l6) . 

From the current density j(x,y) the magnetic induc- 
tion B(.T,y) = (B x ,B y ) = (dA/dy,—dA/dx) is obtained 
via the Biot-Savart law or by using Eq. (21) to get 
A(x,y) = Aj + A a . In strip geometry the magnetic field 
lines are simply the contour lines of A(x, y). In the Lon- 
don case inside the superconductor these field lines co- 
incide with the contour lines of the current density since 
A = —[ioA 2 j, see Figs. 3 and 4. 

Figures 3 and 4 show that the current density exhibits 
a sharp finite peak in the corners; for the depicted case 
a = b = 40A this enhancement is j(a, b)/j(a, 0) = 3.7 for 
H > 0, and j(a,b)/j(a,0) = 5.2 for / > 0. For small A, 
this enhancement increases as A~ 1/3 . Clearly, this cur- 
rent peak favors the nucleation of vortex loops Q| at the 
corners and thus reduces the penetration field. It may 
also enhance the nonlinear Meissner effect 

Figure 5 shows the magnetic moment M(a, 6, A) of 
London strips with various aspect ratios b/a as func- 
tion of A and normalized to the ideal Meissner moment 
M(a,b, 0), Eq. (6), and to the known M(a, oo, A) of the 
infinite slab. One has 
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with the limiting cases 1 — A/aforA^;a and a 2 / (3 A 2 ) 
for A>(i. For thin strips with b <C a and b < A we find 



the limits (see also the Ohmic strip and disk in Ref. |18|) 
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A 2 < ab , 
A 2 > ab . 



Note the nonanalytic A dependence of M for small A, 
which can be seen with the curve b/a = in Fig. 5. 

In summary, we found the exact analytical solution for 
the magnetic moment and surface screening currents of 
long strips with rectangular cross section in the ideal- 
screening Meissner state generated by a homogeneous 
magnetic field H and/or transport current /. Accounting 
for a finite London penetration depth A, we present some 
explicit limiting expressions and numerical results which 
show a high and sharp but finite peak of the current den- 
sity j(x, y) along the four corners. This sharp peak favors 
the penetration of magnetic vortices from the corners, in 
form of quarter loops spanning the corner. From the 
known j(x, y) the exact shape and growth of these loops 
can be obtained in principle, and thus both the micro- 
scopic Bean-Livingston barrier |^,|3| and the macroscopic 
geometric barrier Q-Q] , as well as the thermally activated 
penetration of vortex loops from the corners can be in- 
vestigated in detail. 

Finally, this pronounced current peak is expected to 
enhance the nonlinear Meissner effect predicted in e?-wave 
superconductors |q] and to explain its dependence on the 
sharpness of the specimen corners as observed in Ref. [Q| . 
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FIG. 1. The Meissner surface currents J(x,b) and J(a,y), 
Eqs. (4,5), in strips with various aspect ratios b/a, see inset. 
The dashed line gives the thin strip limit, Eq. (12). 




FIG. 2. Magnetic moment of Meissner strips versus the as- 
pect ratio b/a, Eq. (6) (solid line). The solid line with dots 
shows the fit (8), and the dashed line the limits (9) and (10). 




FIG. 3. The current density j(x,y) along a London strip 
with square cross section (a = b) and London depth 
A/a = 0.025 in a perpendicular magnetic field H. A quar- 
ter of the cross section is shown. Note the sharp but finite 
peak in the corner. The inset shows the magnetic field lines. 




FIG. 4. The current density j(x, y) along the London strip 
of Fig. 3 (a — b = 40A) but with transport current / and no 
applied field, H = 0. 
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FIG. 5. Reduced magnetic moment of London strips 
M* = M(a, b, X)/M(a, b, 0) for various aspect ratios b/a. Top: 
Plotted versus a scaled London depth (\/a)Wl + a/b) the 
curves M*(a, b, A) almost collapse. Bottom: M*{a, b, A) re- 
ferred to the longitudinal limit M*(a, oo, A), Eq. (26). 
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